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Abstract 



subject to unified boundary conditions from F-theory GUTs. We find that very sim- 
ple assumptions about modulus dominance SUSY breaking in F-theory unification 
lead to a predictive set of boundary conditions, consistent with all phenomenological 
constraints. The second lightest scalar Higgs H 2 can get a mass m# 2 ~ 125 GeV and 
has properties similar to the SM Higgs. On the other hand the lightest scalar Hi, 
with a dominant singlet component, would have barely escaped detection at LEP 
and could be observable at LHC as a peak in Hi — > 77 at around 100 GeV. The LSP 
is mostly singlino and is consistent with WMAP constraints due to coannihilation 
with the lightest stau, whose mass is in the range 100 — 250 GeV. Such light staus 
may lead to very characteristic signatures at LHC and be directly searched at linear 
colliders. In these models tan/3 is large, of order 50, still the branching ratio for 
B s — > fi + fi~ is consistent with the LHCb bounds and in many cases is also even 
smaller than the SM prediction. Gluinos and squarks have masses in the 2 — 3 TeV 
region and may be accessible at the LHC at 14 TeV. No large enhancement of the 
H2 —> 77 rate over that of the SM Higgs is expected. 



1 Introduction 



The starting into operation of LHC is already testing many avenues beyond the 
Standard Model (SM). In particular, the discovery of a boson with mass around 
125 GeV [U|2] and properties compatible with those of the SM Higgs is significantly 
constraining many of these ideas beyond the SM. In this regard one may argue 
that such value for a Higgs mass goes in the direction of low energy SUSY, since 
supersymmetric models predict a lightest Higgs with mass < 130 GeV. On the 
other hand, the observed mass is close to the maximum expected in low energy SUSY 
theories, implying a certain degree of fine-tuning in the SUSY-breaking parameters 
which must be relatively large. This is also consistent with the no observation as 
yet of any supersymmetric particle at LHC. 

The simplest testing ground for low energy SUSY is the Minimal Supersymmetric 
Standard Model (MSSM) which does not involve any new particle beyond the SUSY 
partners required by supersymmetry. Still the MSSM has an unattractive ingredient 
in its bilinear Higgs superpotential term, the /x-term. Although supersymmetric, 
this mass term has to be (for no good reason) of order of the SUSY-breaking soft 
terms to get consistent electro-weak (EW) symmetry breaking and low energy SUSY 
spectrum. Perhaps the most economical solution to this problem is the scale invariant 
Next to Minimal Supersymmetric Standard Model (NMSSM) in which a singlet 
S is added to the MSSM spectrum and the /x-term is replaced in the superpotential 
by new couplings, 

Wnmssm = W Yuk + XSH u H d + ^S 3 . (1.1) 

There are no mass parameters in the superpotential and the role of the /x-parameter 
is now played by \{S) which upon minimization of the scalar potential gets naturally 
of the same order than the SUSY-breaking soft parameters. 

The SUSY-breaking soft terms involving the singlet S and the Higgs chiral fields 
have the general form 

y s S of t = m 2 Hu \H u \ 2 + m 2 Hd \H d \ 2 + m 2 s \S\ 2 + (\A x SH u H d + ^A K S 3 + /i.e.) 

(1.2) 

As in the case of the MSSM, the most general NMSSM model has plenty of free 
parameters. On the other hand, in the presence of some underlying unification 
structure at a high energy scale, one expects the number of parameters to be reduced 
to a few. In the case of the MSSM, the constrained MSSM (CMSSM) has universal 
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parameters m, M, A, ji and B, where m is the universal scalar mass, M the universal 
gaugino mass, A the universal trilinear parameter of the standard Yukawa couplings 
and B the universal bilinear coupling, all of them defined at the unification scale. 
Indeed, there are models based on String Theory that lead to such universal selection 
of parameters or extensions of it (for instance with non- universal Higgs masses), see 
e.g. [5] and references therein. Similarly, a constrained version of the NMSSM is 
usually defined in terms of the five universal parameters 

m , M , A, k , A = A x = A K . (1.3) 

In practice, however, one usually takes as free parameters m, M, A, tan /3 and A 
(plus the sign of the effective //-term). The values of and k are fixed upon 
minimization and hence ms in general does not unify with the rest of the scalars of 
the theory. The theory is therefore not constrained in the same sense as it is in the 
CMSSM [6l[7]. One may argue that the singlet may be a bit special and is perhaps 
not surprising that mg is not unified with the rest of the scalar masses. But then 
it would be inconsistent to unify the trilinear A-terms. In particular, A\ and A K 
should be unrelated to the Yukawa trilinear coupling A. Thus, the least one can 
say is that the partially constrained versions of the NMSSM considered up to now 
are slightly inconsistent, unless one allows A\ and A K as free parameters, with the 
resulting reduction of predictivity. This is one of the main issues that we address in 
this paper, namely we try to understand and constraint the structure of the NMSSM 
parameters at a more fundamental level. 

The NMSSM has received recently a lot of attention after the evidence and 
subsequent discovery of the 125 GeV boson at LHC [8H22]. There are two main 
reasons for that. On the one hand, in a general NMSSM model the mass of the 
Higgs particle receives extra contributions from the \SH u H d superpotential term, 

m\ ~ m 2 z cos 2 2/3 + A Vsin 2 2/3 + 6 {ml), (1.4) 

where v = 174 GeV is the Higgs vev, mz is the mass of the Z boson and ^(m^) 
denote the loop corrections to the Higgs mass. In the MSSM these loop corrections 
account for the increasing of the Higgs mass from 90 to 125 GeV, requiring large 
soft terms, stop mixing and fine-tuning. In the NMSSM, however, the second term 
in eq. (11.41) gives an additional contribution to the mass for relatively large A (> 0.5) 
and small tan/3. This allows to get a fairly heavy Higgs boson while reducing the 
fine-tuning. The second reason for this recent interest on the NMSSM is the fact that 
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for some regions of the parameter space one may get an enhanced Higgs decay rate 
to two photons, as suggested by the ATLAS and CMS results as of July 2012 [H[2]. 

As we have already mentioned above, our investigation is however not led by 
these two interests but rather by the attempt to define a fundamentally-motivated 
constrained NMSSM and to check it against the present experimental data, includ- 
ing a 125 GeV Higgs. Concretely, as in the CMSSM, a unified gauge symmetry 
like SU(5) naturally induces universal gaugino masses and also unifies many of the 
sfermion masses. Further assumptions about the origin of SUSY-breaking may lead 
to an increased degree of unification of the SUSY-breaking soft terms. That is for 
instance the case of modulus dominance SUSY-breaking in F-theory SU(5) GUTs 
(for reviews see ref. [23]) that we consider here. In refs. [24"]|25] boundary conditions 
of the general form0 

m 2 = ^|M| 2 , A = --M , (1.5) 

were phenomenologically analysed in the context of the MSSM. These conditions 
appear naturally in F-theory SU (5) schemes in which one assumes that the auxiliary 
field of the local Kahler modulus T is the dominant source of SUSY-breaking, see [M] • 
In the above references it was found that these boundary conditions are consistent 
with all low energy constraints, including a Higgs field with mass around 125 GeV 
and appropriate dark matter relic density. The scheme is very predictive, implying 
tan/5 ~ 41 and a relatively heavy spectrum with M ~ 1.4 TeV, leading to squarks 
and gluinos of order 3 TeV. Such heavy spectrum is required after imposing the 
recent bounds on the branching ratio BK(B S — > /i + /U~) from LHCb [26] and CMS 
[27]. Actually, the most recent results for this decay [28] 129] corner very much the 
parameter space of this very constrained MSSM model. 

On the light of the above NMSSM discussion, it is thus natural to explore whether 
the F-theory SU (5) unification idea exploited in [21] for the MSSM may be extended 
to the case of the scale invariant NMSSM. In this paper we find that the boundary 
conditions ( II. 5p are indeed consistent with all the current phenomenological con- 
straints, including a 125 GeV Higgs boson. A large value of tan/3 ~ 50 is again 
selected and a very small A parameter is required. Small values of A naively corre- 
spond to an effective MSSM limit. Still, the presence of the singlet S leads to quite 
different physics as compared to the MSSM case. First, the lightest supersymmetric 
particle is mostly singlino and correct dark matter abundance is obtained due to 
coannihilation with the lightest stau, which is the NLSP, with masses in the range 
1 Gauge fluxes may slightly distort these boundary condition as we will see later [24] . 
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100 — 250 GeV. Secondly, in addition to the SM-like Higgs with mass around 125 
GeV, there is a lighter neutral scalar with mass around 100 GeV which would have 
barely escaped detection at LEP and that should be detectable at LHC. Thirdly, due 
to an interference effect, the branching ratio BR(I? S — > is easily consistent 

with the recent LHCb data and may be even smaller than the SM prediction. Finally, 
the squark/gluino spectrum is typically lighter than in the MSSM, with masses as 
low as 2 TeV, easily accessible at LHC at 13 TeV. All in all, we find that these are 
highly constrained and predictive scenarios which pass all the current experimental 
constraints with a scale of SUSY-breaking that can be as low asM~ 850 GeV. 

Another particularly attractive aspect of our theoretical approach is the addi- 
tional information obtained about the singlet sector. In the context of modulus 
dominated SUSY-breaking within local F-theory SU(5) unification, the singlet sec- 
tor is less determined. Still, in the simplest scenarios we expect additional constraints 
that, to first approximation, are of the form 

A x ~ —M , A K ~ m| ~ . (1.6) 

Due to the approximate nature of these constraints, we do not impose them in 
our analysis. However and remarkably, we find that almost all the region of the 
parameter space that the current low energy data selects is also consistent with the 
relationships (jl.fip at the unification scale. This is particularly remarkable for the 
linear relationship A\ ~ — M, see figure [10] in section HI Moreover, the geometric 
structure of the F-theory unification predicts small values for the couplings A and 
k, in agreement with the phenomenological requirements. Thus, at least at the 
semi- quantitative level, combining eqs. ( Il.5p and (jl.6p . the present NMSSM scheme 
essentially depends on only three free parameters, namely M, A and k, that are 
reduced (approximately) to two parameters M and A once we impose correct EW 
symmetry breaking. Of course, this is only semi-quantitative since eqs. (11. 6p are 
only approximate and eqs. (II. 5p and (II. 6p may have small corrections from gauge 
fluxes in the extra dimensions. However, it explains the very constrained structure 
of possible benchmarks in this scheme, as displayed in table [1] in section |5j 

Finally, let us mention that in this constrained version of the NMSSM we do not 
find an enhanced Higgs decay rate to 77 compared to the SM, but rather a slight 
reduction. Whereas our numerical analysis takes into account all the known two 
loop radiative corrections to the Higgs sector of the NMSSM [30], in this region of 
the parameter space these corrections are very large. Since some of the two loop 
corrections to the Higgs sector in the NMSSM have not yet been computed and are 



4 



thus not taken into account by our numerical analysis, we expect relatively large 
uncertainties in our numerical estimations for the masses and decay rates of the 
Higgs sector. Nevertheless, a strong enhancement in this rate, as suggested by the 
preliminary ATLAS and CMS recent results PQE], appears to be clearly disfavoured 
in this scenario. 

The structure of this paper is as follows. In section [2] we describe the F-theory 
unification structure in the context of NMSSM models with Kahler modulus domi- 
nation. In section |3] we run the SUSY-breaking soft parameters from the unification 
scale down to the EW scale according to the renormalization group equations and 
perform a computerized scan over the complete parameter space of these models. 
We also comment on the main phenomenological constraints that shape the allowed 
regions of the parameter space. The main constraints arise from the Higgs sector, 
the dark matter relic density and rare decays such as B s — > fi + n~. In section H] we 
analyze the consistency of the regions of the parameter space that pass all the current 
phenomenological constraints, with the F-theory unification structure discussed in 
section |2J In section Owe give some benchmarks and discuss their possible signatures 
at LHC and other future experiments. Finally, we end with some last comments in 
section [6j All the above analysis takes into account the known two loop radiative 
corrections to the Higgs sector of the NMSSM. In appendix |A] we repeat the analysis 
of the Higgs sector but keeping only account of the one loop radiative corrections to 
this sector, showing still qualitatively good agreement with the results in the main 
part of the text. 

2 The NMSSM and F-theory unification 

F-theory SU(5) unification [31] has been recently the subject of intense research both 
from the phenomenological and string compactification points of view [23, 32J . These 
theories provide an ultraviolet completion within string theory to more traditional 
SU(5) GUTs, while involving new mechanisms to address some of the problems of 
field theory GUTs. The symmetry breaking from SU(5) down to the SM is induced 
by the generic presence of hypercharge fluxes in the compact dimensions. The pres- 
ence of these fluxes can also give rise to doublet-triplet splitting for appropriate flux 
choices. Furthermore the observed departure of Yukawa couplings for D-quarks and 
charged leptons from unification may also be accounted for [33]. The reader not 
interested in the formal details of this class of model may safely jump to section [3J 
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F-theory may be considered as a non-perturbative generalization of type IIB 
orientifold compactifications. Ten- dimensional type IIB string theory contains a 
complex scalar dilaton field r = + iCo, where <ft controls the perturbative loop 
expansion and Co is the Ramond-Ramond scalar. The theory is invariant under 
a SL(2,Z) symmetry generated by the transformations r — > 1/r and r — > t + 
i. F-theory provides a geometrization of this symmetry by adding two (auxiliary) 
extra dimensions with T 2 toroidal geometry and identifying the complex structure 
of this T 2 with the complex dilaton r. The resulting geometric construction is 12- 
dimensional and N = 1 4d vacua can be obtained by compactifying the theory 
on a complex Calabi-Yau (CY) 4-fold X±. Such 4-fold is required to be an eliptic 
fibration over a 6-dimensional base B3, so that locally X4 ~ T 2 x 7? 3 . Codimension-4 
singularities of the fibration correspond to 7-branes wrapping 4-cycles of the base 
B 3 . In F-theory 5{7(5) GUTs one set of such 7-branes, wrapping a 4-cycle 5, yield 
an 5{7(5) gauge symmetry. Moreover, for suitable topologies of B 3 , it is possible 
decouple the local dynamics associated to the 5{7(5) branes living on the 4-fold S 
from the global aspects of the B 3 compact space. 

Chiral matter appears at the complex 1-dimensional pairwise intersections of 7- 
branes, corresponding to an enhanced degree of the singularity (see figure [1]). In 
F-theory language the locus of the intersection is usually called matter curve. In 
minimal 5{7(5) GUTs the gauge symmetry is locally enhanced to SU(Q) or 50(10) 
at the matter curves. Indeed, recalling the adjoint branchings 

SU(Q) -> SU(5) x U(l) (2.1) 

35^ 24 + 1 + [5i + c.c] 

50(10) ->• SU(5) x U(l)' (2.2) 

45^ 24 + 1 + [10 4 + c.c] 

we observe that in the matter curve associated to a 5 (or a 5) representation of 
5{7(5) the gauge symmetry is enhanced to 5{7(6), whereas in the one related to 
a 10 of 5{7(5) the gauge symmetry is enhanced to 50(10). In order to get chiral 
fermions there must be non- vanishing fluxes along the {7(1) and {7(1)' symmetries. 

In addition to the above matter curves, a third matter curve with an enhanced 
£{7(6)' symmetry is required to obtain Higgs 5-plets. Yukawa couplings appear at 
the point-like intersection of the Higgs matter curve with the fermion matter curves, 
as illustrated in figure [TJ At the intersection point the symmetry is further enhanced 
to 50(12) in the case of 10 x 5 x 5h down Yukawa couplings and to E§ in the case 
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Figure 1: General structure of a local F-theory SU(5) GUT. The GUT group lives 
on 7-branes whose 4 extra dimensions beyond Minkowski wrap a 4-cycle S inside 
a complex 3-fold B3, on which the 6 extra dimensions of String Theory are com- 
pactified. Gauge bosons live in the bulk of S whereas quarks, leptons, and Higgsses 
are localized in complex curves inside S. These matter curves (denoted as 10, 5, 
5h and 5h in the figure) correspond to the intersection of the 7-branes wrapping 
S with other U(l) 7-branes, as depicted in the figure of the right hand side. There 
is one matter curve for each SU(5) representation. At the intersection of matter 
curves with Higgs curves 5h and 5h, Yukawa couplings develop (figure taken from 
ref. [32]). 

of 10 x 10 x 5h up Yukawa couplings. 

In order to make contact with SM physics, the SU (5) gauge symmetry must be 
broken down to £77(3) x £77(2) x £7(1). In these constructions there are no massless 
adjoints to make that breaking and discrete Wilson lines are also not available. 
Nevertheless, one can still achieve such breaking by means of an additional flux Fy 
along the hypercharge generator in £77(5), leading to the same symmetry breaking 
effect as an adjoint Higgssing. Interestingly enough, this hypercharge flux is also a 
source for doublet-triplet splitting of the Higgs multiplets 5# + 5#. 

Besides the above points of £0(12) and E e enhancement, F-theory £77(5) GUTs 
may also contain points at which the gauge symmetry is enhanced to £77(7). These 
correspond to 5 x 5 x 1 intersections between two 5-plet curves and a singlet curve, 
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Figure 2: General structure of a local F-theory SU (5) GUT with NMSSM couplings. 
The matter curves where the Higgs multiplets live intersect with a transverse curve 
where the singlet S lives. At the triple intersection point the gauge symmetry is 
enhanced to SU(7). 

as can be derived from the adjoint branching 

SU{7) SU(5) x U{1) 2 (2.3) 
48 ^ 24 + l ( o,o) + l(o,o) + [5(o,-i) + l(-i,i) + 5 (1)0) + c.c. ] 

The bosons associated to the extra U(l)s are in general anomalous and become 
massive in the usual wayE 

The simplest models of this class contain one extra singlet S coupling to the 
Higgs multiplets and at low energies are equivalent to the scale invariant NMSSM. 
In what follows we consider that specific setup. The structure of the superpotential 
is thus given by 

Wnmssm = W Yuk + \SH u H d + ^S 3 . (2.4) 

2 F-theory models with hypercharge flux GUT group breaking and SU(5) singlets charged under 
U(l) Peccei-Quinn like gauge symmetries have been argued to require the presence of exotics in 
the spectrum [5Hl37j . Here we assume models with no exotics or with very massive ones. Other 
constructions, such as those where the singlets are not charged under the extra U(l)s (for instance, 
if they come from closed and/or open string moduli in the microscopic theory) also fit well within 
this context. 
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with W-Vuk the standard superpotential Yukawa couplings. To get this structure the 
matter curves associated to the two Higgs 5-plets H u and must intersect with a 
transverse singlet S curve, see figure [2j At the point of intersection the symmetry is 
enhanced from 577(6) to 577(7). Identifying the singlet in the brackets of eq. ( 12. 3p 
with S we see that indeed it may couple to the two Higgs multiplets while respecting 
the U(l) symmetries. 

In an F-theory construction as above there is a Kahler modulus T& whose real 
part tb describes the size of the 6-dimensional manifold B%. In addition, the volume 
of the 4-cycle S where the SU(5) degrees of freedom live is controlled by a local 
modulus T with real part t. Finally, the £7(l)s live in other 4-cycles different than 
S whose volumes are controlled by Kahler moduli that we will denote collectively 
as Ts, with real part t$- In a local setting one assumes t, t$ <C tb, so that to leading 
approximation we can safely ignore the details of B 3 in order to study the dynamics 
of the SU(5) branesfl One can describe this structure by a Swiss cheese type Kahler 
potential of the simplified form [39] 

K = -flog (if - t 3/2 - tf) (2.5) 

where we are working in Planck units. Following [10] (see also [5]), we can make 
use of scaling arguments in order to extract the leading moduli dependence of the 
Kahler metrics of the matter fields. For a matter field X living on a divisor Sx the 
Kahler metric admits an expansion for tx <C tb of the form 

K x = ^— (2.6) 

tb 

where £x is the so-called modular weight of the field with respect to the local modulus 
tx- For matter fields localized on a matter curve the mentioned scaling arguments 
yield £x = 1/2. In our case we have respectively for the SU(5) matter multiplets 
and the singlet S 

+1/2 .1/2 

#5,io = — , K s = (2.7) 
In particular, note that the Kahler metric of S is independent of t. This fact will be 
relevant for the structure of soft terms discussed below. 

The values of tb and t are related to known physical quantities. In particular in 
the perturbative regime one obtains (see e.g. [5]) 

tt = Z U) 9l /3 , t = aa 1 ~ 24 (2.8) 



3 The existence of a manifold admitting such a local limit is however not always guaranteed 
(see e.g. [38]). 
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where «g ~ 1/24 is the gauge coupling at the unification scale, M s = (tV) l l 2 is the 
string scale and g s = (e^) the string coupling constant. For M s ~ 10 16 GeV one has 

1 /3 

tft ~ 115 g s ^> t for g s ^> 0.01, consistently with our requirement t -C £&. 
2.1 The couplings A and k 

Due to the above hierarchy of volumes, physical Yukawa couplings (namely, those 
that correspond to normalized kinetic terms) may differ substantially from the holo- 
morphic Yukawas appearing in the superpotential. Physical Yukawa couplings /i p h ys 
are related to the holomorphic ones ho by the standard N = 1 supergravity expres- 
sion 

= hfiK^K^V (2.9) 

where Ki are the (diagonal) Kahler metrics of the fields involved in the coupling 
and we do not sum over i, j and k. In particular, for the NMSSM specific Yukawa 
couplings in eq. (12 .4p we obtain 

1/2 

A phy s = Aot- 1/2 t 5 1/4 ~ A ^, (2.10) 

,-3/4 
K phys — K % 

where we have made use of eqs. (I2.5P and (12.71) . For ts — t ~ 24 we thus get a 
suppression factor in relating the holomorphic couplings Ao and Ko to the physical 
ones, with A p h ys , K p h ys — 9 x 10~ 2 Ao, ^o- We will show in section |3] that appropriate 
EW symmetry breaking in these models together with the current experimental 
constraints require small physical couplings A p h ys and K p h ys - This is indeed consistent 
with the above suppression coming from the Kahler metrics, as typically |Aq|, |«o| — 
1. On the contrary, large values for these couplings as required e.g. in refs. [8HT0], 
[T2Tnfil[r8H2T)| 22J do not seem viable within our context. 

We have already mentioned in the previous subsection the possible origin of the 
cubic A coupling in eq. (12.41) from the intersection of three matter curves on a point 
of enhanced SU(7) symmetry. If that is the case, the cubic k coupling in eq. (12. 4p 
should arise from instanton corrections that violate the (typically anomalous) U(l)s, 
under which 5* is charged (see e.g. [H]). In particular, the absence of a quadratic 
S* 2 coupling suggests the presence of a Z 3 gauged symmetry remnant of a gauged 
U(l)x symmetry [42]. One thus expects k to be a small parameter 

4 Cubic self-couplings could also arise if S were in fact a linear combination of three different 
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2.2 T-modulus dominance and unified soft-terms 



A natural source of supersymmetry breaking in type IIB string theory compactifica- 
tions, or more generally in F-theory, is the presence of certain classes of closed string 
antisymmetric fluxes (G4 fluxes in F-theory). As described in |43j . from the point of 
view of the 4-dimensional effective supergravity, the supersymmetry breaking fluxes 
are encoded in non-zero vevs of the F-auxiliary fields of the Kahler moduli. In this 
work we assume a hierarchy of vevs F tb 3> F t ^> F ts as it will lead to a very con- 
strained set of soft-terms. Such hierarchies of auxiliary fields easily arise in large 
volume models (see [39]). Thus, following [24J and making use of eqs. (I2.5P and (12 .7p 



we can compute explicitly the structure of soft terms at the unification scale for the 
MSSM sector of the theory, yielding results as in [24J, namely 



M = j , (2.11) 
2 \M\ 2 A 3 



\M\ 



2 



A = -y(3-pfr) • 

where M are universal gaugino masses, A are trilinear couplings (which appear 
multiplied by the SM Yukawa couplings in the Lagrangian) and m#, mg and mio are 
universal masses for the scalar fields in the three SM matter curves. The parameter 
Ph corresponds to a correction describing the effect of gauge fluxes on the Higgs 
matter curve, see [21] for details. This parameter should be small, of order pu ~ 
a l Q 2 ~ 0.2 or smaller. 

In addition to the above MSSM-like soft terms, there are also soft terms that 
involve the singlet S. The structure of such soft terms is however more subtle 
since, unlike the SM fields and gauge bosons, 5* is not localized on the 4-cycle S 
and therefore may be subject to extra sources of supersymmetry breaking. The 
statements that one can make for the soft terms involving S are thus more model 
dependent. In the simplest case, with no other sources of supersymmetry breaking 
for S other than F tb and F t , making use of the Kahler metric for S, eq. (12 .7p . yields 
at the unification scale 

A x = -M(l-p H ), A K = m| = 0, (2.12) 



singlets within an extended gauge symmetry £7 or Eg, in F-theory. Even in this case a small 
coupling k is expected due to the geometric suppression discussed above. 
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where A\ and A K are the coefficients of the trilinear scalar couplings X(SH u Hd) and 
kS 3 in the Lagrangian .3 

In the next section we analyze the phenomenological implications of the bound- 
ary conditions (12. lip that depend on the two parameters M and pu, and we leave 
A, k, A\, A K and mj as free parameters. However, as described in detail in section 
HI the regions of the parameter space that are consistent with correct EW symmetry 
breaking, a Higgs mass in the vicinity of 125 GeV and other phenomenological con- 
straints tend to select solutions that are also consistent with the constraints (12. 12ft ! 
This is remarkable since in no place such second stronger unification conditions were 
imposed. 



3 Higgs masses, dark matter and other constraints 

In order to analyze the phenomenological prospects of the NMSSM models with F- 
theory unification boundary conditions discussed above, in this section we run the 
soft terms eqs. (I2.1ip from the unification scale down to the EW scale according to 
the renormalization group equations, and impose radiative EW symmetry breaking 
in the standard way as well as the main phenomenological constraints. The strongest 
phenomenological constraints that shape the structure of our solutions come from the 
Higgs sector, the dark matter relic density and the branching ratio BK(B S — > p + p~), 
as we discuss below. 



3.1 Scanning over modulus dominated NMSSM vacua 

Without imposing the extra conditions ( I2.12p . NMSSM models with F-theory unifi- 
cation boundary conditions consist of seven parameters, namely the values of M, pn, 
A, k, A\, A K and ms at the GUT scale. One combination of these parameters is how- 
ever fixed by requiring the correct pole mass for the Z boson, mz = 91.187 GeV [H]. 
Equivalently, we can take the six independent parameters to be given by the values 
of M, p H , A x and A K at the unification scale, A at the supersymmetry breaking 
scale and tan /3 at the scale mz- We perform a large scan over this parameter space 
by means of a modified version of the computer program NMSSMTools [15] v3.2.0, 
where we have implemented the boundary conditions eqs. ( 12. lip at the unification 

5 By a slight abuse of notation, we make use of the same letter to denote the multiplet (appearing 
in the superpotential) and its scalar component (appearing in the Lagragian). 
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scale. We take the latest combined result from Tevatron for the top quark pole mass, 
m t = 173.2 GeV [46]. 

The output of NMSSMTools consists of the supersymmetric mass spectrum com- 
puted at the scale Q 2 = m^m i3 , with q 3 and u 3 the third generation squarks; 
as well as the main couplings and reduced cross sections at the supersymmetry 
breaking scale, taking into account the known two-loop radiative corrections to the 
NMSSM Higgs sector [30]. Moreover, it imposes the main experimental constraints 
from LEP, Tevatron, LHC and WMAP, that we have updated with the most re- 
cent CMS results for the observed 95% confidence level upper limit on the reduced 
cross sections for the Higgs decays H — > t + t~ [47], H — > W + W~ — > t^t'vv [48J, 
H ZZ ->■ £+£+£-£- [49\ and VH Vbb |5DJ, based on 17 fb" 1 of integrated lumi- 
nosity at y/s = 7 - 8 TeV, the ATLAS and CMS recent results for H -> 77 [51]l52]. 
based on near 11 fb _1 of integrated luminosity at =7 — 8 TeV, and the AT- 
LAS result for the MSSM decay H + — > t + v [53], based on 4.6 fb _1 of integrated 
luminosity at yfs = 7 TeV. We have also included the recent LHCb result mea- 
suring BR(B S — > = 3.2^2 x 1CT 9 at 95% confidence level [29]. Regarding 
the relic density of neutralino dark matter, NMSSMTools calls the computer program 
micrOMEGAs [5l] v2.2 in order to estimate it and compares with the amount of cold 
dark matter observed by the WMAP satellite, 0.1008 < Vth 2 < 0.1232 at the 2a 
confidence level [55] . We do not impose a priori any constraint on the muon anoma- 
lous magnetic moment (g — 2) M but will comment on the numerical results for this 
quantity below. 

Despite the large amount of constraints, we find solutions with correct EW sym- 
metry breaking that pass all the current theoretical and phenomenological con- 
straints and have a moderately low scale of supersymmetry breaking, M > 850 
GeV. All these vacua have large tan/3 and small A, typically with tan/3 > 46 and 
A < 0.01 (see figure |3] below and figure [8] in section 144]) . The input parameters for 
a representative sample of benchmark points are shown in table [TJ where we also 
include k and rrig UT . The resulting values of the NMSSM parameters at the 
SUSY scale are indicated in table [2j Having such large values of tan f3 introduces 
important uncertainties in the computation of the masses for the Higgs sector, as 
threshold corrections to down Yukawa couplings become large in this regime (see 
e.g. [SUEZ] for similar considerations in the context of the MSSM). Moreover, in 
this region of the parameter space also two-loop radiative corrections to the neutral 
scalar and pseudoscalar masses are large and particularly sensitive to the value of 
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Point 


M 


tan j3 


PH 


x gut 


K GUT 


A<3UT 


J^GUT 


rnGUT 


Pi 


858.2 


47.5 


0.035 


5.2- 10~ 3 


4.1- lO" 4 


-895.4 


-319.2 


105.9 


P2 


888.0 


49.0 


0.0092 


6.5- 10~ 3 


3.8- 10- 4 


-922.9 


-195.6 


50.0 


Pz 


981.4 


49.5 


0.021 


6.7-10- 3 


4.1- 10- 4 


-1015.7 


-262.2 


82.9 


Pi 


1009.9 


52.9 


0.195 


6.4- 10- 3 


3.7- 10- 4 


-1045.7 


-252.7 


78.5 


P 5 


1036.0 


51.5 


0.135 


6.3- 10~ 3 


3.9- 10- 4 


-1070.3 


-303.1 


99.9 


Pe 


1180.5 


53.3 


0.030 


5.5-10- 3 


2.5- 10- 4 


-1204.2 


-198.2 


50.6 


P7 


1218.5 


56.8 


0.30 


5.7- 10" 3 


2.8- 10" 4 


-1257.9 


-281.6 


91.9 


Ps 


1271.0 


55.1 


0.182 


5.1- lO- 3 


1.9- 10- 4 


-1290.0 


-160.5 


22.7 



Table 1: F-theory unification/modulus dominance boundary conditions for a repre- 
sentative of benchmark points satisfying all the experimental constraints. Masses 
are given in GeV. 
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3.8- 10- 4 
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49.9 
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P 3 


6.1- lO" 3 


4.1- 10~ 4 
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82.8 


1401 


Pi 


5.7- lO" 3 


3.7- 10- 4 


-98.9 


-252.6 


78.4 


1447 


P5 


5.7- lO" 3 


4.0- 10- 4 


-107.3 


-303 


99.9 


1478 


P 6 


4.9- 10- 3 


2.5- 10~ 4 


-91.0 


-198.2 


50.5 


1662 


P7 


4.9- lO" 3 


2.8- 10- 4 


-107.3 


-281.5 


91.9 


1723 


Ps 


4.5- 10- 3 


1.9- 10- 4 


-86.0 


-160.5 


22.5 


1780 



Table 2: Boundary conditions at SUSY scale. Masses are given in GeV. 

m t . Whereas most of these radiative corrections have been computed [30j and are 
actually taken into account by NMSSMTools, some of them remain yet unknown. The 
reader should hence bear in mind that some of our numerical estimations (mostly 
those concerning the masses of the neutral scalar and pseudoscalar sectors) may 
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contain relatively large uncertainties^. Nevertheless, we find that the qualitative 
conclusions of our analysis are solid, and do not change if the analysis is done for 
instance at the one loop level (see appendix |A|) . 

On the other hand, having small values of A is consistent with the F-theory 
considerations discussed in section 12.11 and drives these solutions near the effective 
MSSM limit of the NMSSM. In spite of the apparent similitude with the MSSM, 
there are however important differences, coming from the fact that for such small 
values of A the lightest scalar and pseudoscalar Higgsses and the lightest neutralino 
are largely dominated by their singlet and singlino components, respectively. This 
has important consequences for the neutralino relic density and the Higgs couplings, 
as we discuss below. 

Whereas we have scanned for arbitrary signs and values of A\, A K and A, we have 
only found phenomenologically viable vacua for the sign choice A\ < 0, A K < and 
A > 0. We will discuss further this point in section 14.21 In what follows we focus 
in that region. We have represented in figure |3] the values of A and k at the scale 
of supersymmetry breaking for the set of NMSSM vacua considered in this paper. 
As we observe, for A > we recover positive values of k and generally one order 
of magnitude smaller, k < A/10. This is a consequence of the minimization condi- 
tions of the scalar potential and is consistent with what happens in the constrained 
NMSSM 0. 



3.2 The Higgs sector at two loops 

In the NMSSM the Higgs sector consists of three neutral scalars, H iy i = 1,2,3, 
two neutral pseudoscalars, Ak, k = 1,2, and one charged scalar, H + . Scanning 
over modulus dominated NMSSM vacua as described in the previous subsection, 
reveals that the boundary conditions ( 12. lip do not allow for > 120 GeV in 
any region of the parameter space. The only possibility is hence that H 2 plays the 
role of the SM-like Higgs boson, with m# 2 > m^, and H\ is fairly decoupled from 
the SM sector such that it escaped detection at LEP. This leads naturally to small 
values of A since in that case if i is mostly singlet (see figure E]) and it couples very 
weakly to the SM fields. Thus, one of the main predictions of NMSSM models with 
F-theory unified boundary conditions (12.111) is the presence of an extra light Higgs 

6 A similar statement is obtained from comparison with the results obtained by means of the 
SARAH-Spheno [58] NMSSM implementation presented in [59 that includes the two-loop radiative 
corrections of [30] . 



15 



A 




Figure 3: Distribution of NMSSM vacua over the plane (A, k) with F-theory boundary 
conditions at the unification scale given by eqs. (|2.1ip and unconstrained values of A\, A K 
and A. Dark blue and red points both pass all the current experimental and theoretical 
constraints. Dark blue points, however, have a deficiency of neutralino relic density and 
therefore require some additional source of dark matter. 

that is fairly decoupled from the SM fields. Such possibility has also been recently 
considered in [2D]. 

To be more precise, we represent in figure @] the distribution of masses (m Hl ,m H2 ) 
for the two lightest Higgs bosons. The strongest constraints in this plane come from 
the LHC upper limits on the reduced cross section of the Higgs decay H — > 77 
[51], 152"] and from the analogous LEP results for H — )■ bb [50] ■ The former puts a 
lower bound m# 2 > 122 GeV on the mass of the SM-like Higgs, whereas the latter 
constrains vacua where the lightest Higgs H\ couples too strongly to the SM fermions. 
Since in these models the scalar singlet component is distributed among H\ and H2, 
and large masses m# 2 require a non-negligible singlet component for H 2 , the cross 
section a bbZ (e + e~ ->• H X Z) increases with the mass of H 2 and therefore LEP bounds 
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Figure 4: Distribution of the light Higgs masses (mH 1 ,mH 2 ) for NMSSM vacua with F- 
theory unification boundary conditions (|2.1ip and unconstrained values of A\, and A. 
Dark blue and red points both pass all the current experimental and theoretical constraints. 
Dark blue points, however, have a deficiency of neutralino relic density and therefore 
require some additional source of dark matter. 



effectively constrain m# 2 from above, see figure HI This upper bound on m# 2 is 
in particular (slightly) stronger than the one derived from the recent LHC results. 
Vacua that are consistent with both LHC and LEP bounds thus have approximately 
rriH 1 = 100 ± 15 GeV and m# 2 = 124 ± 2 GeV, with the signal of Hi fitting in the 
2a excess observed at LEP [60J. 

It is also interesting to compare the predicted reduced cross sections of the SM- 
like Higgs boson H 2 with the recently observed Higgs signal at the LHC, particularly 
in the H — y 77 channel, as it is starting to being measured with increasing precision. 
In figure[5]we represent the reduced signal cross section in the gg — > if 2 — > 77 channel 

°sm(99 -> H) 



17 




Figure 5: Distribution of the the reduced signal cross section in the H 2 — > 77 channel, 
{99 ~~ H2) versus the mass for the SM-like Higgs H 2 in NMSSM vacua with F-theory 
unification boundary conditions (|2.11|) and unconstrained values of A\, A K and A. Dark 
blue and red points both pass all the current experimental and theoretical constraints. 
Dark blue points, however, have a deficiency of neutralino relic density and therefore 
require some additional source of dark matter. 



against the mass of H 2 . As we have already mentioned, the singlet component of 
H 2 is not negligible for the above allowed range of masses. This leads to a mild 
suppression of the couplings of the SM-like Higgs H 2 to the other SM fields, and 
in particular to the top quark and the W bosons that dominate the one loop SM 
contribution to a^i^gg — > H). Moreover, the stau is not light enough to enhance the 
di-photon production by running in the loops. Hence, as it can be observed in figure 
there is no enhancement in a 77 (g g — > H 2 ) with respect to the SM but rather a 
mild suppression, with R 2 r< {gg — > H 2 ) ~ 0.7 — 0.9. Although this is still in reasonable 
agreement with the latest LHC results [5UE2], an experimental confirmation of a 
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large enhancement in K^{gg — > H 2 ) could disfavour the present class of modelsJ3 
For the same above reasons, similar considerations apply to the other reduced cross 
sections of H 2 , leading also to mild suppressions with respect to the SM. Concretely, 
we observe Rf +W ' (gg H 2 ) ~ R% z {gg -> H 2 ) ~ R T 2 +T ~ (VBF -»■ # 2 ) ~ 0.7 - 0.8. 

We now briefly comment on the remaining part of the Higgs spectrum. Apart 
from Hi and H 2 discussed above, the pseudoscalar A\ is also relatively light in these 
models, with rriA 1 ^ 350 GeV, see figure |6j However, it is highly dominated by its 
singlet component (its doublet composition is insignificant, of the order of ~ 10~ 8 %) 
and therefore very decoupled from the SM fermions. In this limit the pseudoscalar 
mass can be written as m\ w — 3kv s A k . On the contrary, H%, A 2 and H + have 
much larger masses and are nearly degenerated, with m# 3 ~ m^ 2 ~ itl h + > 1 TeV. 
They are therefore much harder to detect at the LHC. 

Finally, we may wonder about the robustness of the above Higgs spectrum. As 
we have already mentioned, in this region of the parameter space of the NMSSM 
there are relatively large uncertainties coming from some of the two loop corrections 
to the NMSSM Higgs sector that have not yet been computed. In this regard and 
for completeness, in appendix |A] we perform the same analysis of this subsection 
but keeping only one loop radiative corrections to the Higgs sector into account. 
Such analysis reveals that the above qualitative results for the Higgs sector also hold 
at the one loop level, although the range of masses is broadened considerably and 
small enhancements to the di-photon rate with respect to the SM appear also to be 
possible. 

3.3 Neutralino dark matter 

The lightest neutralino, Xii can be a viable dark matter candidate in the NMSSM 
with interesting phenomenological properties [01]. It differs from the MSSM in 
that it now contains a singlino component, which alters its couplings to the SM 
particles. The neutralino is also very sensitive to the structure of the Higgs sector, 
since it determines its annihilation cross section in the Early Universe (and thus 
the theoretical predictions for its relic density) and plays an important role in the 
computation of the scattering cross section of quarks. 

Nevertheless, as noted in appendix \Kl mild enhancements in the di-photon rate with respect to 
the SM are allowed by keeping the computation of the Higgs masses at the one loop level. In this 
regard it would be interesting to understand the effect of the complete set of two loop corrections 
to the NMSSM Higgs sector on the di-photon rate. 



19 



m H2 (GeV) 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 

50 100 150 200 250 300 350 



m Al (GeV) 

Figure 6: Distribution of the SM-like Higgs and the lightest pseudoscalar masses 
(mH 2 > m A\ ) for NMSSM vacua with F-theory unification boundary conditions (|2.1ip and 
unconstrained values of A\, A K and A. Dark blue and red points both pass all the current 
experimental and theoretical constraints. Dark blue points, however, have a deficiency of 
neutralino relic density and therefore require some additional source of dark matter. 

In our construction, the neutralino turns out to be an almost pure singlino in all 
the allowed regions of the parameter space. This is due to the very small values of 
n, which are at least one order of magnitude smaller than A as we already showed 
in figure [3J This leads to a hierarchical structure in the neutralino mass matrix, 
kv s <C M 12 < Xv s , which implies that the gaugino and Higgsino components of 
the lightest neutralino are almost negligible (to less than approximately a 0.1%). 
This also implies that the neutralino mass can be small in these scenarios without 
violating any experimental limit. To a good approximation we can write m^o w 2kv s , 
and in principle, neutralinos as light as 50 GeV are possible, but when the recent 
experimental constraints are applied on other observables (mostly the Higgs decays 
discussed above), we are left with the range 100 GeV < m^o < 240 GeV. 
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We have represented in figure [7J the mass of the lightest neutralino m^o against 
the stau mass m.f, for vacua that pass all current experimental constraints except 
that of the relic density observed by WMAP. A singlino-like neutralino has reduced 
couplings to the SM particles and thus generally displays a small annihilation cross- 
section in the Early Universe. As a consequence, the relic density is too large in most 
of the parameter space, and typically exceeds the WMAP range 0.1008 < Qh 2 < 
0.1232 (at the 2a level). The relic density can be lowered through coannihilation 
effects when the stau mass is very close to that of the neutralino. In fact, we observe 
that the allowed regions of the parameter space have m? — m%a ~ 10 GeV. This 
leads to a very interesting prediction of this scenario, namely that the stau mass is 
in the range 110 GeV < m f < 250 GeV. Notice that unlike the scenarios discussed 
in ref. [22], the stau- neutralino mass difference is never small enough to allow for 
long-lived staus. 

Dark matter can be detected directly through its scattering off nuclei in a de- 
tector. The interaction of neutralinos with quarks has contributions from s-channel 
squark exchange and t-channel Higgs exchange diagrams. Since the neutralino is a 
pure singlino state, its coupling to squarks is completely negligible. A pure singlino 
only couples to the singlet part of the Higgs, but this coupling C§ ss = — 2y/2n van- 
ishes in the limit k — >■ and is tiny in our vacua. This results in a extremely small 
neutralino-nucleon scattering cross-section. In particular, the theoretical predictions 
for the spin-independent contribution is below ~ 10 -13 pb, several orders of 

X\ v 

magnitude below the predicted sensitivity of the projected 1 ton scale detectors. The 
same happens with the spin-dependent component. 

3.4 Constraints on rare decays and the muon anomalous 
magnetic moment 

Rare decays constitute excellent probes for new physics beyond the SM. In particu- 
lar, the effect of flavour-changing neutral currents in 6-physics signals is extremely 
interesting in supersymmetric models, since it can be sizeable at large tan /3, and 
rather sensitive to the Higgs sector. For example, supersymmetric contributions to 
the branching ratios of the rare processes B s — > fi + fi~ or b — > sj can easily exceed 
the experimental measurements of these quantities and generally lead to stringent 
constraints on the parameter space. This is particularly important in the NMSSM 
due to the presence of new scalar and pseudoscalar Higgsses which induce new con- 
tributions to these observables. 
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Figure 7: Distribution of neutralino and stau masses (m^o,mf) for NMSSM vacua with 
F-theory unification boundary conditions (|2. 1 1 [) and unconstrained values of A\, A K and 
A, for points that pass all current experimental constraints except that of the relic density, 
that is encode in the legend of colors. 



Let us start by addressing the branching ratio BR(i? s — > fi + fi~). This observ- 
able can be written in terms of the Wilson coefficients which appear the effective 
Hamiltonian that describes the transition b — > s as follows, 



BR(B S 



fx H 



Ami 



m 



B 



2m, 



m 



B 



(3.2) 



In the SM calculation of this quantity [62J only Ca is relevant, since Cs and Cp are 
suppressed by the small Yukawas. However, in supersymmetric theories, there are 
penguin contributions involving the neutral scalar and pseudoscalar Higgs bosons to 
Cs and Cp which can be sizeable, both in the MSSM [631 and the NMSSM 



In the MSSM one finds BR(5 S — > fi + fi ) oc (tan 6 (3/m A ), and therefore in the large 
tan {3 regime supersymmetric contributions easily exceed the recent LHCb measure- 
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ment BR(B S — > fi + fi~)= 3.2^ 2 x 10~ 9 [29], especially for light pseudoscalars. In 
the NMSSM the Wilson coefficient Cp receives contributions from both pseudoscalar 
Higgses, but it is only their doublet component that contributes. 

As we have already commented in section I3T2"} the viable points resulting from our 
scan display a relatively light pseudoscalar, rriA 1 < 350 GeV, which is a pure singlet. 
Therefore, it does not contribute to the Wilson coefficient Cp, and only the heavy 
pseudoscalar A2 has to be taken into account. The latter is rather heavy (in our scan 
its mass varies in the range from 1 to 1.4 TeV) and as a consequence the resulting 
value of Cp is also small. In fact, it is of the same order as the SM contribution 
(2mn/m 2 B ) Ca but with opposite sign. This implies that not only BK(B S — > 
can be small, in agreement with the experimental data, but even a cancellation 
between Cp and Ca that leads to a smaller value than the SM prediction is possible 
(see [57J for a similar effect in the context of the MSSM). This last possibility, in 
fact, happens in parts of our parameter space for M > 1200 GeV, where values as 
low as BR(i? s — > w 1.8 x 10 -9 can be obtained. As we will see in the next 

section (see figure E]), for not too large values of tan/3 this actually favours a light 
spectrum. 

We will now briefly turn our attention to the branching ratio BR(6 — > 57). In the 
past years this has been one of the strongest constraints on supersymmetric models, 
mainly because of the additional contributions from loops of charged Higgs bosons. 
There are also specific contributions from the extended Higgs and neutralino sectors 
of the NMSSM, although these start at the two loop level [64"tl65]. Whereas imposing 
the experimental result BR(5 — > S7) = (3.52 ± 0.23 ± 0.09) x 10~ 4 generally leads to 
constraints on the NMSSM parameter space [57], our scan reveals that these bounds 
are now superseded by those on BR(B S — > fi + fi~) discussed above. Therefore, the 
current experimental bounds on BR(6 — > 57) do not have a large impact in the space 
of parameters of these models. 

Let us finally address the supersymmetric contribution to the muon anoma- 
lous magnetic moment, a^ USY . The observed discrepancy between the experimen- 
tal value [68] and the Standard Model predictions using e + e~ data, favours posi- 
tive contributions from new physics. These can be constrained to be in the range 
10.1 x 10~ 10 < a s ^ USY < 42.1 x 10~ 10 at the 2 a confidence level [59] where theoret- 
ical and experimental errors are combined in quadrature (see also refs. (TUIITT] that 
lead to a similar range). However, if tau data is used this discrepancy is smaller 
2.9 x 10 -10 < a^ USY < 36.1 x 10~ 10 [71]. The theoretical predictions for this observ- 
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Table 3: Low energy observables and reduced cross-sections for different Higgs signals 
calculated for the set of benchmark points. 

able generally increase for light supersymmetric masses and large tan/3. Although 
in this scenario the masses of the first and second family sleptons are relatively high, 
we are in the regime of large tan/3. The predicted a^ USY turns out to be within the 

3 a range of the e + e~ data and well within the 2 a range of the discrepancy predicted 
by tau data. 

To illustrate this discussion, we display the specific values of these low energy 
observables for our choice of representative benchmark points in table [31 

4 Consistency with F-theory unification 

In the previous section we have discussed the main phenomenological constraints that 
shape the parameter space of NMSSM vacua with F-theory unification boundary 
conditions eqs. (12. lip . We now describe the overall structure of the allowed regions 
and further comment on their consistency with an underlying F-theory unification 
structure, including also the extra conditions f)2.12p . 

Indeed, as we have already mentioned, leaving the singlet sector unconstrained, 
NMSSM vacua with F-theory unification boundary conditions are completely spec- 
ified by six parameters that encode the relevant information of the local F-theory 
background. In our scan these are given by the values of M, p H , A\ and A K at the 
unification scale, A at the supersymmetry breaking scale and tan at the scale mz- 
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However, in specific F-theory GUTs the singlet sector will also satisfy a set of bound- 
ary conditions at the unification scale which, under the simplifying assumptions of 
section [21 are given by eqs. (I2.12p . Hence, the number of independent parameters in 
specific models can be reduced from six to three (for instance, the values of M, pu 
and A at the unification scale), as we discuss below. 



4.1 M, p H and tan (3 

Let us first address the effect of the phenomenological constraints on the region of the 
parameter space spanned by M, p H and tan We have represented in figure [8] the 
distribution of phenomenologically viable NMSSM vacua with F-theory unification 
boundary conditions (12. lip , over the relevant region of the plane (M, tan 0). As we 
have advanced in previous sections, vacua satisfying all the current phenomenological 
constraints sit in the region of large tan with tan j3 > 46. In this region the 
contribution of the heaviest pseudoscalar A 2 to the branching ratio BR(B S — > 
is such that it leads to important cancellations with the SM contribution, giving rise 
to comparable or even reduced values with respect to the SM (see section I3.4p . 
For regions of the parameter space where the above cancellation is maximal the 
branching ratio BR(B S — > fi + p~) in particular becomes too small, below the recent 
95 % confidence level lower limit measured by the LHCb collaboration [29] . For 
moderately large values of tan this disfavours a band in the plane (M, tan 0) with 
approximately 1200 GeV < M < 1800 GeV, see figure Nevertheless, the reader 
should bear in mind that all points in this band have BR(B S — > p + pr) > 1.8 x 10 -9 , 
and therefore many of them are actually within the theoretical uncertainties from 
QCD [63]. 

Regarding the scale of supersymmetry breaking, we observe in the same figure 
that the lower bound on M compatible with all the current experimental constraints 
is M ~ 850 GeV. Below that scale, the cross section for the SM-like Higgs decay 
i?2 — > t + t~ becomes larger than the most recent experimental upper limits set by 
the CMS and ATLAS collaborations [4THT2] . On the other hand, since the LSP 
is a singlino-like neutralino, its mass m^o ~ 2ks < 240 GeV is independent of M 
and vacua with the correct neutralino relic density are possible for arbitrarily large 
values of M, contrary to what occurs in the modulus dominated MSSM [24, 25J. 
Nevertheless, the amount of fine-tuning required in order to have a light stau with 
mass in the region of co-annihilation increases considerably for large values of Mjf 



^Moreover, for large values of M the value of tan (3 becomes too large for the computation of 
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Figure 8: Distribution of NMSSM vacua over the plane (M, tan /3) with F-theory bound- 
ary conditions at the unification scale given by eqs. (|2.1ip and unconstrained values of A\, 
A K and A. Points pass all the current experimental and theoretical constraints. Dark blue 
points, however, have a deficiency of neutralino relic density and therefore require some 
additional source of dark matter. 

Similarly, we have represented in figure M the distribution of phenomenologi- 
cally viable vacua on the plane (pn, tan j3). In this case, the LHCb upper limit on 
BR(B S — > forces pu to take small values, pu < 0.4, with most of the points 

below 0.2. Remarkably, this is consistent with our treatment of the effect of mag- 
netic fluxes on the Kahler metrics for the Higgs as small perturbations, see section 

13 

Having looked at the effect of the phenomeno logical constraints on M, pn and 
tan (3, let us now move on to the discussion of the singlet sector. 

down-type Yukawa couplings to be reliable. 
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tan (3 BR(B S n + fjr) < 2.0 x 10" 9 (LHCb '12) 



Figure 9: Distribution of NMSSM vacua over the plane (ph, tan f3) with F-theory bound- 
ary conditions at the unification scale given by eqs. (|2.1ip and unconstrained values of A\, 
A K and A. Points pass all the current experimental and theoretical constraints. Dark blue 
points, however, have a deficiency of neutralino relic density and therefore require some 
additional source of dark matter. The upper yellow and pink regions in this plot are only 
indicative: whereas all points inside this regions are excluded because the branching ratio 
BR(i? s — > p + p~), mn% and/or Qh 2 , there are also points outside this regions that are also 
excluded by these observables. 

4.2 The singlet sector: A\, A K and A 

As already mentioned in section [2J the singlet S is not localized in the local 4-cycle S 
and therefore its interactions depend on the details of the background along the bulk 
of the compactification, making this sector very mo del- dependent. Nevertheless, in 
simple models where the only source of supersymmetry breaking for S is given by the 
F-terms F t and F tb we expect A\, A K and to satisfy a set of boundary conditions 
at the unification scale which, to first approximation, are given by eqs. f )2.12p . 
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It is highly remarkable to see how the above type of relations are also strongly 
favoured by the current experimental input. Indeed, although we have scanned over 
unconstrained values and sign choices of A\, A K and A, requiring a SM-like Higgs in 
the window 122 — 128 GeV and a neutralino relic density below the WMAP bound 
flh 2 < 0.1232 uniquely selects the sign choice A\ < 0, A K < and A > 0, consis- 
tently with the F-theoretic boundary conditions for the singlet sector, eqs. ( I2.12p . 
Moreover, it is possible to check that the current experimental constraints imply 
also a strong correlation between the values of A\ and M at the unification scale, 
similar to that of eqs. (12.121) . Indeed, in figure [10] we have represented the distribu- 
tion of NMSSM vacua on the plane (M, A\) with boundary conditions ( 12.111) and 
unconstrained values of A\, A K and A. The linear correlation between A\ and M 
for vacua that pass all the current phenomenological constraints is clearly visible 
in this figure. Finally, in all these vacua A K and ms take relatively small values at 
the unification scale, with —0.015 < (rns/ m 5,io) 2 ^ 0.035 and (A K /A) < 0.3 (see 
tabled] in section [5]), also in fairly good agreement with the F-theoretic expectation 
A K ~ m s — in eqs. (12.121) . 

In order to make a more quantitative estimation of the compatibility between 
the phenomenological constraints and the F-theoretic relations for the singlet sector, 
we can consider a slightly more general Kahler metric for 5* than that of eq. (12 .7p . 
depending also on the Kahler modulus of S, 

Ks = t -^— (4.1) 

Making use of this metric, the boundary conditions ( 12 .121) generalize to 

A x = -M(2-Z t -p H ) (4.2) 
A K = -3M(1 - &) 
™! = |M| 2 (l-6) 

where the simple F-theory unification conditions for the singlet sector, eqs. fl 2 . 1 2 j) . 
are recovered for £ t = 1. Making use of these equations we can estimate the modular 
weight £ t in three different ways (one from each equation) for each vacuum in our 
scan with boundary conditions (12.111) and unconstrained values of A\, A K and A. 
If the vacuum is consistent with the F-theoretic boundary conditions, the three 
estimations of £t should agree with each other and with £ t ~ 1. We have represented 
in figure [11] (left) the maximum value of the three estimations against the mass of 
the SM-like Higgs H 2 for vacua that pass all the current experimental constraints. 



28 



M (GeV) 



r 



■ ■ ■ 



1800 — 




1000 



1600 — 



1400 — 



1200 — 




800 — 



-1800 



-1600 



-1400 



-1200 



-1000 



A x (GeV) 



Figure 10: Distribution of NMSSM vacua over the plane (M, A\) with F-theory boundary- 
conditions at the unification scale given by ecjs. (12. lip and unconstrained values of A K 
and A. Points pass all the current experimental and theoretical constraints. Dark blue 
points, however, have a deficiency of neutralino relic density and therefore require some 
additional source of dark matter. For convenience, in this plot we have not imposed the 
lower bound BR(£ S -> n + fi~) > 2.0 x 10~ 9 . 

We observe in this figure that the latter strongly favour the F-theory value £t = 
1. Similarly, in figure dH (right) we have represented the maximum percentage of 
mismatch between the three different estimations of £ t , showing that for a large 
fraction of the phenomenoiogically viable vacua the mismatch between the three 
eqs. (14. 2 p is lower than 10%. 

Hence, whereas in the scan over modulus dominated vacua that we have per- 
formed in section 13.11 we have taken six independent parameters, it turns out that 
many of the vacua that satisfy all the current experimental bounds can be actually 
described to a large extend in terms of just three parameters, namely M, pu and 
A, satisfying the boundary conditions (12. lip together with the additional boundary 
conditions (I2.12p for the singlet sector. 
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Figure 11: (Left): Modular weight £t estimated from eqs. (j4.2[) against the mass of the 
SM-like Higgs, for NMSSM vacua with F-theory boundary conditions at the unification 
scale given by eqs. (|2.1ip and unconstrained values of A\, A K and A. (Right): Percentage 
of mismatch between the three eqs. (|4.2|) against M. For convenience, in this plot we have 
not imposed the lower bound BR(.B S — > fi + fi~) > 2.0 x 1CP 9 . In both figures all the points 
pass all the current experimental and theoretical constraints. Dark blue points, however, 
have a deficiency of neutralino relic density and therefore require some additional source 
of dark matter. 



5 Supersymmetric spectrum and signatures 

The model is extremely predictive. After imposing all the experimental constraints 
on the parameter space we are left with a very characteristic supersymmetric spec- 
trum. In tableH] we indicate the masses of the various supersymmetric particles and 
Higgs bosons for a series of representative benchmark points. 

Regarding the Higgs sector, this model provides a very interesting scenario, with 
two light scalar Higgses. In this sense it seems similar to the proposal of ref. [20] 
where the lightest Higgses had masses of 98 and 125 GeV. Notice however that 
the solution that we find belongs to a completely different region of the NMSSM 
parameter space. The most important differences are the large value of tan (3 in our 
scenario and the very small values of both A and k. Also /i e // > 1 TeV, thus being 
much larger than in ref. [20]. This leads to a different phenomenology (as already 
emphasized in section l3~2| most notably, in our case we predict no enhancement 
of Higgs decay H 2 — > 77). Notice also that this scenario does not contain exotic 
channels for Higgs decay of the two lightest CP-even states, since the pseudoscalar 
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is heavier than Hi and H 2 . In any case, we can be certain that this scenario can be 
tested through the branching fractions of these two lightest states. In particular, the 
lightest Higgs could lead to a peak in the Hi — > 77 channel around 100 GeV that 
would be observable at the LHC. Regarding the lightest pesudoscalar, it is singlet 
like and its production rate is extremely suppressed. The heavier pseudoscalar and 
the third CP-even Higgs state are doublet like, but very heavy and thus very difficult 
to produce. 

The coloured section is in the 1.5 to 3 TeV range. Squarks are lighter than 
the gluino and, interestingly, the lightest stop can be as light as ~ 1200 GeV 
in the regions with smaller gaugino masses. As already emphasized the next-to- 
lightest supersymmetric particle is the lightest stau, whose mass difference with the 
neutralino is of the order of 10 GeV for the whole range of viable gaugino masses. The 
nature of the neutralino and chargino sector is easily understood from the resulting 
hierarchy kv s <C Mi < M2 < Xv s . As already stated in section l3~3| the lightest 
neutralino is pure singlino. The second neutralino state is bino-like and relatively 
light {fn^o 300 — 500 GeV). The third neutralino state and the lightest chargino 
are wino-like and with a mass in the range m-o ^± ~ 600 — 1000 GeV. Finally 

A3 1X1 

the heavier neutralino states and the heavier chargino are Higgsino-like states with 
masses approximately equal to /i e // = Xv s . 

This kind of spectrum is very similar to what is obtained in the Constrained 
NMSSM [6] and both share the same search strategies. In particular, the branching 
ratios of the lightest chargino and the second and third neutralino states into the 
stau NLSP are sizable, thus potentially leading to significant rates of tau-rich final 
states [TUlEl]. Let us be more specific about this. Left-handed squarks can decay 
into the lightest chargino with a branching fraction of approximately 65% and into 
the third neutralino with 33%, whereas right-handed squarks decay mostly into the 
second neutralino with a branching fraction of approximately 99%. Regarding the 
stops, the lightest stop undergoes the following decays ti — > (60%), ii — > x®t 
(12%), ti x^t (25%). So either if they are directly produced of obtained in 
gluino decays, the resulting production of xt an d X23 * s copious. The second and 
third neutralino states decay in turn into the lightest stau with a branching ratio 
of approximately 100%. Obviously the stau NLSP can only decay into the lightest 
neutralino, producing another tau in the final state which is softer than the previous 
one (as the mass- difference of the stau and neutralino is approximately 10 GeV) 
thus leading to ^2,3 ~~ > tit ~~ ^ T+T ~Xv The lightest chargino decays into a stau and 
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Point 


9 


Qr,l 


*1,2 


&1,2 




n,2 




xi 


m-Hi 


m Ai 


m H + 


Pi 


1921 


1758 
1827 


1263 
1558 


1481 
1583 


684 
827 


230 
719 


213 

367; 696 
1238; 1243 


696 
1244 


103 

122.4 

1016 


321 
1016 


1019 


P 2 


1983 


1814 
1886 


1302 
1601 


1521 
1626 


708 

855 


175 
735 


164 

381; 721 
1274; 1279 


721 
1279 


98.1 

123.4 

1036 


220 
1036 


1040 


Ps 


2716 


1989 
2069 


1434 
1749 


1673 
1778 


782 
944 


199 
807 


189 

423; 800 
1394; 1400 


800 
1400 


96.9 

124.8 

1131 


273 
1131 


1134 


Pi 


2236 


2042 
2125 


1499 
1802 


1718 
1827 


804 
971 


197 
831 


186 

436; 824 
1440; 1444 


824 
1445 


97.4 

124.3 

1095 


266 
1094 


1098 


P 5 


2289 


2091 
2175 


1527 
1841 


1762 
1868 


825 
996 


216 
851 


205 

447; 845 
1471; 1475 


845 
1476 


97.4 

124.7 

1148 


306 
1148 


1151 


Pe 


2585 


2358 
2455 


1728 
2064 


1986 
2095 


939 
1133 


178 
955 


167 

513; 967 
1653; 1657 


967 
1657 


98.5 

124.4 

1274 


223 
1274 


1277 


Pi 


2663 


2428 
2528 


1809 
2134 


2046 
2161 


970 
1169 


204 
990 


193 

530; 999 
1712; 1716 


999 
1716 


93.9 

123.4 

1227 


287 
1227 


1230 


Ps 


2769 


2525 
2629 


1862 
2207 


2127 
2238 


1011 
1219 


164 
1023 


153 

554; 1043 
1770; 1774 


1043 
1774 


97.9 

123.7 

1330 


192 
1329 


1332 



Table 4: Supersymmetric spectrum and Higgs masses for the set of benchmark 
points. All the masses are given in GeV. 

the corresponding neutrino thus giving rise to only one final state tau xt ~~ >* v t^\ — > 
u T TXi- The signal expected from this kind of scenarios is therefore the presence of 
multitau signals, originated from the two chains of cascade decays, associated to the 
emission of hard central jets and missing energy [73J. 

Notice finally that the upper bound set by BR(5 S — > fi + n~) implies that the 
whole spectrum is lighter than approximately 3 TeV. This is well within the reach 
of LHC at 14 TeV for searches involving multijets plus missing energy. 
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6 Discussion 



The most elegant solution to the /i-problem in the MSSM is perhaps its extension 
to the scale invariant NMSSM model. In the MSSM, the existence of unifying 
underlying symmetries suggest the existence of unified boundary conditions leading 
to simple structures like that of the constrained MSSM, the CMSSM, with only a few 
free parameters, i.e. m, M, A, /i and B. In contrast, unlike the MSSM case, there is 
some ambiguity in the definition of what a constrained version of the NMSSM, the 
CNMSSM, could be. 

In this paper we have analyzed in detail a constrained version of the NMSSM 
model with boundary conditions obtained from the assumption of modulus domi- 
nated SUSY breaking in F-theory SU(5) unification models. Such theories provide 
for an ultraviolet completion for more traditional SU(5) GUTs and provide for new 
solutions for problems like doublet-triplet splitting and D-quark to lepton mass ra- 
tios. One obtains a simple and very predictive structure of SUSY breaking soft terms 
summarized in eqs.f l2.ll"]) . Additional approximate boundary conditions eq. fll.6p may 
also be derived for soft terms involving the NMSSM singlet S. Furthermore, the ge- 
ometrical structure of F-theory unification implies that the NMSSM couplings A and 
k are small with A, k < 0.1. In spite of all these very constrained parameter values, 
we find that the obtained very constrained NMSSM model is consistent with correct 
EW symmetry breaking, a 125 GeV Higgs and appropriate relic dark matter. The 
model passes also B-decay constraints as well as LEP and LHC limits. 

Independently of its underlying string theory motivation, the structure of soft 
terms f)2.1ip and f |1.6p provides for a definition of a constrained NMSSM model 
consistent with all presently available data. The resulting NMSSM model is very 
predictive and the low-energy spectrum very constrained, as it is visible from table 
(JH). The second lightest Higgs scalar H 2 can get a mass m# 2 ~ 125 GeV whereas 
the lightest scalar Hi, with a dominant singlet component, barely escaped detection 
at LEP and could be observable at LHC as a peak in Hi — > 77 at around 100 GeV. 
The LSP is mostly singlino and may provide for the correct relict density thanks to 
coannihilation with the lightest stau, which is the NLSP with a mass in the range 
100-250 GeV. Such stau NLSP leads to signals at LHC involving multi-tau events 
with missing energy. In this model the value of tan f3 is very large, tan /3 ~ 50, still 
the branching ratio for B s — > fi + fi~ is below recent LHCb and CMS bounds and in 
many cases is smaller than the SM prediction, due to an interference effect. Gluinos 
and squarks have masses in the 2 — 3 TeV range with a stop with a mass as low as 
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1.2 TeV, all within reach of LHC(14). On the other hand no large enhancement of 
the if 2 — > 77 over that of the SM is expected. It is exciting to think that, if indeed 
the signatures above are observed at LHC, it will be evidence not only for SUSY 
but for an underlying F-theory unified structure. 
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itih 2 (GeV) 




m Hl (GeV) 

Figure 12: Distribution of the light Higgs masses (mH 1 ,mH 2 ) for NMSSM vacua with 
F-theory unification boundary conditions (|2.1ip and unconstrained values of A\, A K and 
A, computed at one loop. 

A The Higgs sector at one loop 

For completeness, in this appendix we repeat the analysis of section 13.21 taking only 
into account one loop radiative corrections, as well as the two loop contributions 
to the bottom and top Yukawa couplings to leading logarithmic approximation. 
Namely, in this appendix we do not consider two loop radiative corrections to the 
Higgs sector. As we summarize below, the results however agree qualitatively with 
the ones in section 13. 2\ showing the robustness of our conclusions. 

In the models that we consider in this paper two loop radiative corrections tend 
to lower Higgs masses in as much as 2 - 5 GeV. Thus, at one loop we find in 
general higher masses than those in section 13.21 In spite of this, the mass of the 
lightest Higgs, Hi, is still not large enough to fit the LHC Higgs signal at any 
region of the parameter space, so that the role of the SM-like Higgs boson has 
again to be played by H 2 . In figure [T2] we represent the distribution of masses 
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rriH 2 (GeV) 




Figure 13: Distribution of the the reduced signal cross section in the H2 — > 77 channel, 
Rip (gg — > H2), versus the mass for the SM-like Higgs H2 computed at one loop in NMSSM 
vacua with F-theory unification boundary conditions (I2,ll|) and unconstrained values of 
A x , A K and A. 

(m^m^) f° r the two lightest Higgs bosons computed at one loop. One difference 
with respect to the two loop analysis performed in section I5T21 is that the LHC bound 
on a w+w {gg — > H) now becomes also relevant, excluding solutions in the region 
m Hl > 110 GeV, m H2 > 129 GeV. Moreover, the LEP bound on a blz (e + e- ->• HZ) 
and the LHC bound on a" f ' y {gg — > H) become softer, leading to the range of masses 
50 GeV < m Hl < 120 GeV and 123 GeV < m H , 2 < 129 GeV. The upper limit 
m A 1 ^ 350 GeV on the mass of the lightest pseudoscalar, on the other hand, still 
holds. 

Regarding the reduced cross sections of the SM-like Higgs boson H 2 , we find that 
the contribution from stau loops to the diphoton rate production become relevant 
in some cases, allowing for small enhancements with respect to the SM, as it can 
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be seen in figure [THJ More precisely, we observe from that figure that the reduced 
signal cross section l£P(gg — > H 2 ) lays in the range 0.6 — 1.4. 
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